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£\j ' Abstract 

The long-range correlation energies and the off-diagonal interactions are studied and a general formula for correlation 
energy E c of the 7r electron systems is given, which is beyond the nearest-neighbor electron-electron interactions and includes 
, the off-diagonal interactions. It is found that the effects of the off-diagonal interactions W and X on the correlation energies 
■ are opposite, but the influence of X on the correlation energies is much stronger than that of W on the correlation energies, 
Q\ ' and the correlation energies decrease with increasing the screening effect. 
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T Introduction 

S the previous paper (Ref.[l]), the author deduced a simple expression for the long-range electronic correlation energies E c 
erf the 7r electronic systems where only the nearest-neighbor electron-electron interaction v between two adjacent lattices are 
i-6©*lsidered, and the formula was applied to the polymer polyacetylene (PA) and used to calculate the correlation energies and 
the calculating results are well consistent with those obtained by other theoretical calculations including ab initio method [2]. 
In order to obtain more general correlation energies formula which includes the more long-range Coulomb interactions (that is 
rli.e diagonal interactions) and also the off-diagonal interactions, the correlation energy formula in the Ref.[l] should be revised. 
£N| As we all have known, besides the diagonal interactions, there are off-diagonal interactions in the 7r conjugated polymers. The 
rrrtarest-neighbor diagonal interactions are the on-site Coulomb interaction U (Hubbard interaction) and the nearest-neighbor 
^^ulomb interaction v. The nearest-neighbor off-diagonal interactions have X and W, which are so-called, respectively, the 
CWjnd-charge interaction and the bond-bond interaction [3]. There have been many studies about the roles of the off-diagonal 
Cwe[ractions[4]-[8]. There are also the studies about the effects of the off-diagonal interactions on the excitonic binding energies 
3if the 7r conjugated polymers [9] [10]. It was also known that the off-diagonal interactions have relationship with the screening 
ifom the 7r electrons in the conjugated polymers[ll]. At the different screening and different conditions, the effects of the 
erff-diagonal interactions on the physical properties of the electron systems may have different influences. For instance, X can 
used to describe the possible superconducting states of the organic polymers or the non-7r electron systems [12], and W is 
y^ed to describe the ferromagnetism in the polymers [13]. Until now, these issues have not been studied and there is no detailed 
Studies about the contribution to the correlation energies from the off-diagonal interactions. Thus the correlation energies from 
the off-diagonal interactions should be studied. So, it is necessary to deduce a formula of the correlation energies where the 
off-diagonal interactions have been included explicitly and in the same time the Coulomb interactions are not limited to the 
nearest-neighbor interaction v. 

Therefore, in this study we mainly consider the general situation where both the long-range electron-electron interactions 
and the off-diagonal interactions are considered at the same time as well as the relation of the screening and the correlation 
energies. 



2 Theoretical formation 

As the procedure of deducting the desired formula is similar to that used in the Ref.[l], so here some of the deduction procedure 
are omitted and the detailed are mainly referred to the Ref . [1] . We directly start to our discussion from the following correlation 
energy (see the expression (10) in the Ref.[l]) 

E * = \ E/ d 3 rd 3 r' P<7 (r)v(f~ ?)p% a , (r, ?) = 1(1- I HF ) (1) 
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where v(f — r') = , the Coulomb interaction between two charges at r and r' . In calculation of conjugated polymers, it is 
parameterized by the Ohno potential. Here I equals 

I = Y.J J d 3 rd 3 r'p a (r)v(r- P)p a ,(P)g aa ,{r, P) (2) 

a a 1 

and 

J2 I I d 3 rd 3 r'p a (r)v(r-P)p a/ (P)g^(r,P) (3) 



'HF 

aa' 



Here the pair-distribution g(f,r')(= J2aa' g a a' (r, r'j) satisfies the sum rule: / d 3 r' p(r')[g(f,r') — 1] = —1. Using the relation 
Pa{r) = Pij^ii^jir) and the sum rule, (2) reads 

^EE^ / / d^r'tfi^^vir-P^UPOMr') ~ ^ £ £>« / d 3 rci*(r)v(r)<j>j(r) (4) 

erer' ij,fcZ " " a ij 

The two-fold integration in the first term in (4) is a many-center integral which describes the Coulomb interaction between 
two 7r electrons lying, respectively, at the positions r with spin a and r' with spin a'; it is denoted generally as V^m- As 
for Vij t ki, generally it is hard for us to analytically calculate this many-centered integral easily and simply. In the CNDO 
(complete neglect of differential overlap) approximation [14], Vi^ki reduces to the Coulomb integrals Vn.kk which is a two- 
center integral. In the situation of the conjugated polymers, Vu t kk is further parameterized into the Ohno potential, that is, 
Va^kk = Vi t k = U/y/l + (Jfi — fk)/ao) 2 , where U is the on-site electron-electron interaction[15]. If the overlap integrals are 
explicitly considered, the two-electron interaction integral Vij.ki in (4) may be approximated as USijSki/ 'yl + {{fij — %i)/co) 2 5 
where Sij(= J d 3 r<j)*(r)(j>j(r)) is the overlap integral, and ffj means the position vector between the site i and the site j'[16]. 
When the overlap integrals are omitted, Sij becomes and then Vij.ki goes back to the above expression Vu^k- After 
rearrangement, we obtain 

^ i^k ^ k^l ^ i^j.k^l 

1 = 2^^ VuMP ^' k+ 2EE^^' + 2^ £ VwPiiPu 

aa' i,k aa' i,kl aa' ij,kl 

+ ^££^ + ^££«, (5) 

a i a ij 

where the parameter A has been introduced again by the replacement of e 2 by Ae 2 in v(r — r'). And Vi and Jy are 

Vi = - J d 3 r<P*{r)v{r)<P t {r) = -e 2 J d 3 r^^- (6) 



Jij =~ J dM*(r)v{r)^{r). (7) 

Here the integral Vi is the Coulomb interaction of a positive charge (+e) at the position r' = with the tt electron at the 
position r near the site i. This positive charge is in fact just the correlation hole moving together with the 7r electron. The 7r 
electron and the positive charge is at least a lattice distant apart. The integral (i ^ j) is a two-center integral. 

It can been seen from the definitions of them that off-diagonal interactions refer to the overlaps of two tt wave functions lying 
at two different sites [9] Thus they are actually very small compared with the diagonal interactions Vi»,fcfc[17]. Thus, in ordinary 
theoretical studies, only the diagonal interactions Va,kk ar e kept and the off-diagonal interactions are not considered. So only 
the interactions Vu^k in the first term of (5) are kept in the physical study of low-dimensional n electron systems, just like 
that in the CDNO approximation. Similarly, the integral Jik also refers to the overlap between two n electronic wave functions 
respectively lying at two different lattice sites and is a very small quantity compared with Vu y kk, then it can be omitted in some 
actual theoretical calculations. It can be kept in theoretical calculations where more accuracy is required. 

By these considerations, (5) reduces into the form 

A A 

T - 2EE«^ + 2EE^ ( 8 ) 

aa' i,k a i 

Now we calculate Ihf- Because the procedure of calculating Ihf is given in the Ref.[l], here is not given out and readers 
can be referred to that paper. Finally we have 



i^k 

Ihf = ^V^Mp^-plM. 



(9) 



Table 1: The relations of the off-diagonal interactions and the screening factor when U = 6 eV. The unit is eV. 
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Inserting (8) and (9) into (1), the expression of E c is obtained 

E * = ^El^S + X>^) 2 ] (io) 

cr.i a,ik 

This is a more general formula for the electronic correlation energy. The correlation energy can be acquired using this formula 
when the averages of the diagonal charge density and the off-diagonal charge density p a ik (i ^ k) between the sites i and k and 
also Vi t k (i ^ k) are known. When all \Vi\ is assumed to take the nearest-neighbor Coulomb interaction v and Vi.k — Va,kk = v 
with k = i + 1, the above formula (10) reduces to Eq.(15) obtained in the Ref.[l] which is a special case of this formula. 
When all kinds of electron-electron interactions are kept, the correlation energy formula is given by 

i^k 

a,i aa' i,k 

k^l ^ i^j-k^l ^ i^j 

+ 2EE X »."^« + 2E E WijMPijpH + zYtT,^^ ( n ) 

crcr' i,kl era' ij,kl o ij 

where Xa^i = Vu^i {k ^ I), Wij t ki = (* i>k =/= I). Here the off-diagonal interaction Xa_ki is called the bond-site 

interaction, the off-diagonal interaction Wij^i is called the bond-bond interaction. This is general form of the correlation energy 
which includes all kinds of the electron-electron interactions (the diagonal and the off-diagonal interactions). 



3 Estimation of correlation energies 

For the situation where only the nearest-neighbor interactions are considered, the above (11) can be simplified. For the one- 
dimensional 7r electron conjugated polymers with N carbon atoms such as the bond-alternated chain polyacetylenc (PA), in 
the bond order wave (BOW) phase, the average charge density p^ — 1/2, the bond charge density (also called the bond-order 
matrix) pf i+1 = Pu+i = p + {—l) l 5p- Thus we have the correlation energy formula from (11) 

K = ^[-l + 2(v + W)(p 2 + (Spf) + (X + J)p] (12) 

where v » \X\ > W > 0, v » | J|, and X < 0, J < 0. Note (12) is suitable to an even number conjugated polymer chain with 
single bond and double bond. For a odd number conjugated polymer, there is a more term appearing in (12): — A{v + W)p8p 
and -{X + J)5p. 

The values of v,W and X have been already estimated before for PA, which are shown in Table 1, and the value of J is not 
known. However, we may estimated the value of J through a rational comparison between X and J. According to the definitions 
of the integrals of X(< 0) and J(< 0) in the nearest-neighbor situation, X is J / d 3 rd 3 r'4>*(r)(f)*(r') x v(r — r')^)i(r')4>i+i(r), 
and |J| is about J d 3 rcj)*{r)v{r)(j)i + i(r), where v(f — r') is the Coulomb interaction between two electronic charges, then it 
may be estimated that the maximum of \X\ is about ^- J d 3 r4>*(r)(j>i(f) J d 3 r' '0*(r')0 i+ i(r'), and the maximum of \ J\ is about 
— J d 3 r' <j)*{r')(j)i + i{r'), where a is the lattice constant and the wave functions <pi{r) is normalized. So the values of \X\ and 
|J are approximately comparable with each other, that is, \X\ ~ |J|. Therefore the correlation energies formula for the 
nearest-neighbor interactions reads 

K' = ^-[- V - + 2(v + W){? + {5pf) + 2Xp] (13) 

According to this formula we may estimate the correlation energies under different off-diagonal interactions. In ab innito 
numerical computations for many-electronic systems of atoms, molecules and solid state, the values of correlation energies are 
mixed in exchange-correlation energies E ex instead of a separate quantity, and moreover the off-diagonal interactions can not be 
calculated concretely. Those are shortcomings of the ab inito numerical calculations. Through the formula here we may clearly 
see the effect of the off-diagonal interactions on the correlation energies and may discuss the physical properties such as the 
band structure and electronic structure as well as other physical properties related to elemental excitations. The results of the 
calculations are shown in Table 2. 



Table 2: The correlation energies e c (= E"'/N, N is number of the tt electrons in a conjugated polymer chain) at different 
screening /3 and different off-diagonal interactions X and W when U = 6 eV. The unit is eV. 
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e c (X^0,W = 0) 
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-0.0659 


-0.0597 
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e c (X = 0,W^0) 


-0.1648 
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e c (X^0,W^0) 


-0.1869 


-0.0623 


-0.0539 


-0.0509 



4 Discussion and summary 

The integral V% is an integral of a positive charge (Correlation hole) at the origin with all other electrons in the whole system. 
Although the form of the integral V{ is simple but its calculation still is not direct, and the approximate calculation about it 
is available for some special systems. The simplest is that V{ takes the nearest-neighbor interaction v. The more complicated 
calculation is left in the later study and here we do not discuss its computation. Similarly direct computation of the integral J 
is not considered here and is left in the later study. The site charge density pf 4 (that is, the diagonal elements of charge density) 
and the bond charge density p a ik (i ^ k) (that is, the off-diagonal parts of charge density) can be calculated directly from the 
elliptic integrals for a rigid conjugated polymers [9], and also can be drawn from the results of numerical computation including 
ab initio computations. 

It is all known that the values of the off-diagonal interactions are related to the screening values [3]. Table 1 lists the values 
of the diagonal interaction v and the off-diagonal interactions X and W according to the computing results of the Ref.[3]. 
According to the Ref.[3], the order of the diagonal and the off-diagonal interactions are v > \X\ > W. All interactions (the 
diagonal and the off-diagonal interactions ) are related to the on-site Hubbard interaction which was discussed in the Rcf. [9]. 
When U is given, the values of v, W, and X are determined. In our calculation, [7 = 6. It is seen from the Table 1 that at 
the normal screening case {(3 ~ 1) and the middle screening case {(3 ~ 3), W is much smaller than X, and X is much smaller 
than v. Only with increasing the screening up to the case (3 > 5, both W and X become stronger, they even go over the value 
of v. Anyway at the normal screening and the middle screening, the off-diagonal interactions W and X may be regarded as a 
perturbation relative to v. Interesting is that the off-diagonal interaction X is negative, which may be thought to be related to 
the possibility of the superconductivity of the conjugated polymers[12]. 

It is seen from our calculating results shown in Table 2 that when the screening (3 increases from normal to strong, the 
correlation energies diminish from large to small whatever with the off-diagonal interactions (X, W) and without the off- 
diagonal interaction (X, W). It is seen from Table 2 that the off-diagonal interaction X increases the correlation energies at 
most and the off-diagonal interaction W decreases the correlation energies compared with those without considering the off- 
diagonal interactions X or W. When both X and W are considered in calculations, the correlation energies increase compared 
with those without considering both X and W . It has been known that the bond-site interaction X was used to discuss the 
superconductivity of the polymers [12] [17]. Here the X is shown to have contribution to the correlation energy, which may imply 
important significance of the long-range correlation to the superconductivity of the polymers. 

It has been known that under fixed U, the exciton binding energy decreases with increasing the screening /3[9]. Here under 
fixed U, the correlation energies decrease with increasing the screening (3. About the screening /?, the bond density correlation 
of the 7T electrons, the off-diagonal interactions, the exciton bonding energy, and their relations has been detailed in the Ref.[9]. 
[3 is a cutoff parameter of the Coulomb interaction. The larger the (3, the weaker the Coulomb interaction. Thus the factor 
(3 reflects a kind of screening effect on the n electrons in the conjugated polymers. But the larger the (3, the changes of the 
relative values among v, W, and X, see Table 1. This shows the significance of the screening (3: the larger (3, the stronger the 
bond density correlations. Notice here that we should not confuse the correlation energy with the bond density correlation. 
Therefore, we see that with increasing the screening (3 (> 3), the off-diagonal interaction X increases much, the exciton binding 
energy decreases, the bond density correlation increases, and the correlation energies also increase. It is found from the Table 2 
that the correlation energies decrease rapidly with increasing the screening when (3 is from 1 to 7 when without considering the 
off-diagonal interactions W and X. But when considering only X or considering both X and W, the change of the correlation 
energies with increasing screening {(3 > 3) is not so big, this is because when the screening (3 increases the off-diagonal interaction 
X also increases very much. 

In summary, in this paper the general formula of the correlation energies has been deduced with the off-diagonal interactions 
X and W by the first time. The effect of the off-diagonal interactions X and W on the correlation energies are discussed. It 
is found that the effects of the off-diagonal interactions W and X on the correlation energies are opposite. Given the screening 
(3, the off-diagonal interaction X increase the correlation energies and the off-diagonal interaction W decreases the correlation 
energies. But the influence of X on the correlation energies is much larger than that of W on the correlation energies, and the 
correlation energies decrease with increasing the screening factor (3. 
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